THE DAMPED STRING PROBLEM REVISITED 



FRITZ GESZTESY AND HELGE HOLDEN 



Abstract. We revisit the damped string equation on a compact interval with 
a variety of boundary conditions and derive an infinite sequence of trace for- 
mulas associated with it, employing methods familiar from supersymmetric 
quantum mechanics. We also derive completeness and Riesz basis results (with 
parentheses) for the associated root functions under less smoothness assump- 
tions on the coefficients than usual, using operator theoretic methods (rather 
than detailed eigenvalue and root function asymptotics) only. 



1. Introduction 
We reconsider the damped one-dimensional wave equation 

p(x) 2 u tt (x,t) - u xx (x,t) + a(x)u t (x,i) = 0, (x,t) E (0, 1) x [0, oo), 
u(-,0) = / , «t(-,0) = /i, (1.1) 
it(-,t), t £ K, satisfies certain boundary conditions at x = and x = 1, 
assuming 

< p.p- 1 E L°°([0, l];dx), a £ L°°([0, 1]; dx), a real-valued, (1.2) 

and appropriately choosing fj E L 2 ([0, l];p 2 dx), j = 0, 1. Suitable boundary con- 
ditions for u(-, t), t > 0, in (1.1) at x = and x = 1 studied in this paper are, for 
instance, Dirichlet, Neumann, (anti)periodic boundary conditions, etc. (we refer to 
Section 3 for details). 

We note that the seemingly more general equation with variable speed c(-) > 0, 

p(x) 2 u tt (x,t) - c(x) 2 u xx (x,t) + a{x)u t {x,t) = 0, (x,t) £ (0, 1) x [0, oo), (1.3) 

subordinates to the case described in (1.1) as long as 

< c,^ 1 E L°°{[0,l];dx), (1.4) 

replacing p by p/c 2 and a by a/c 2 in (1.1). 

Our interest in this topic and the principal motivation for writing this paper has 
its origin in a question posed by Steve Cox in September of 2008: He had derived 
the following trace formula (in the case p = 1) 

^L 2 ([o,i]-,dx)(a((-d 2 /dx 2 ) D y 1 ) = dxx(l - x)a(x), (1.5) 

Jo 

in the case where u(-, t), t > 0, in (1.1) satisfies Dirichlet boundary conditions at x = 
0, 1 (in particular, the operator (—d 2 /dx 2 )D in (1-5) denotes the Dirichlet Laplacian 
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in L 2 ([0, 1]; dx)). Steve Cox then posed the question whether there actually exists 
an infinite sequence of such trace formulas in analogy to the well-known sequences 
of trace formulas for completely integrable evolution equations of soliton-type (e.g., 
the Korteweg-de Vries (KdV) equation). This question will indeed be answered 
affirmatively in Theorem 5.4 in the general case where p is nonconstant and for a 
variety of boundary conditions at x = 0, 1. 

We note that the area of damped wave equations remains incredibly active up 
to this day. Since we cannot possibly describe the recent developments in detail in 
this paper, we refer, for instance, to [6], [7], [8], [9], [11], [13], [15], [16], [17], [19], 
[20], [21], [22], [23], [24], [31], [32], [34], [36], [37], [38], [39], [46], [53], [54], [55], [56], 
[57], [58], [59], [60], [66], [74], [75], [76], [77], [78], [79], [86], [87], [88], [90], [91], [92], 
[94], [101], [103], and the references therein, which lead the interested reader into a 
variety of directions. 

The traditional semigroup approach to the damped (abstract) wave equations 
(1.1) (cf., e.g., [33, Sect. VI. 3], [35, Ch. VIII]) consists of rewriting it as a first-order 
system of the type 




lL 2 ([0,l];p 2 dx) 




(1.6) 



where p 2 (-d 2 /dx 2 ) bc > in L 2 ([0, 1]; p 2 dx) is of the type 

p- 2 (-d 2 /dx 2 ) hc =T* c T hc (1.7) 

with 

T bc = ip- 1 (d/dx) bc (1.8) 

in L 2 ([0, 1]; p 2 dx), and the subscript "be" represents appropriate boundary condi- 
tions at x = 0, 1 (Dirichlet, Neumann, (anti)periodic, etc.) to be detailed in Section 
3. 

However, we will not be working with the generator 

_ / h 2 ([0AU P 2 dx)\ 

T -' Q/ " 2 \-p- 2 (-d 2 /dx 2 ) hc -ap- 2 )> (1.9) 

dom(G Tbc , Q/(9 2) = dom(T b * c T bc ) © dom(T bc ), 
in the Hilbert space %T b © L 2 ([0, 1]; p 2 dx) associated with (1.6), where 
^T bc = (dom(T bc ), (•, -)-H Tb )i 

(u,v) Ut ^ = (T hc u,T hc v) L 2 ([01] . p 2 dx) , u,v £ dom(T bc ). 
Instead, we will put the principal focus on the Dirac-type operator 

T,*\ ia I lL 2 (\QA]:n 2 dx) 0\ 



(1.10) 



D hc + B 

' 'be o j P < y o oy (i.n) 

dom(L> bc + B)= dom(L> bc ) = dom(T bc ) © dom(T b * c ) 

in the Hilbert space £ 2 ([0, 1]; p 2 dx) 2 by employing the fact that iGT bc . a / P 2 m %T bc © 
L 2 ([0,l];p 2 dx) is unitarily equivalent to (D bc + B)(I L 2 ([01] . p 2 dx) © P ran(Tbc )) in 
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L 2 ([0, 1]; p 2 dx) 2 (cf. Theorem 2.4 for details), as recently proven in [39] in the 
context of abstract wave equations with a damping term. 

Working with £> bc + B rather than using G Tb . Q / p2 has two distinct advantages: 
First, the unperturbed Dirac operator 




a,c = u, r (i.i2) 



in L 2 ([0, 1]; p 2 dx) 2 is self-adjoint and of a supersymmetric nature (cf. Appendix A 
for details) which permits its spectral analysis in terms of the operator T bc T bc > 
(resp., T bc T bc > 0), and second, the non-self-adjoint term B is represented in terms 
of a simple diagonal 2x2 block operator in L 2 ([0, 1]; p 2 dx) 2 . 

Assuming (1.1) and choosing ( G K\{0} with \(\ sufficiently small, such that 
C £ p(-Dbc + B), we prove in Theorem 5.4 that 

trL2([o,i];p 2 <fa) 2 ( Im [(A>c + £ - ChT 1 ] ) 

= Im [tr i2([04];p2d;E) ((2C + z(a/p 2 ))(T b * c T bc -C 2 /-^(a/p 2 ))- 1 )] 



m=0 



^bc,2mC 2m j (1-13) 



where 

*bc,o - tv L , {%lhp , dx) {{a/p 2 ){T^T hc )- 1 ), (1.14) 
tbc, 2 = -tr L2([0 ,i ]; ^) ((a/p 2 ) (T b * c T bc )- J (a/p 2 ) (T bc T bc ) _1 (a/p 2 ) (T^)" 1 ) 

+ 3tr L2([0 , 1];p2d;c) (( a /p 2 )(T b * c r bc )- 2 ), (1.15) 

etc., 

and explicit boundary condition (be) dependent expressions for ibc,o are listed in 
(5.16)-(5.19). The following infinite sequence of trace formulas, our principal new 
result, is proved in Theorem 5.11: 

^ Im (A J -(£> bc + £?)"'+ 1 ) /-tbc,2n, m = 2n, 

g lA^A. + B)!^) = {o, ro = 2n+l, nGN °- (L16) 

Explicitly, one obtains for m = 0, 1 in (1.16), 

£ ^Ct^F = ~ <bc '° = -^([M^itW^KWc)" 1 ). (1-17) 



x ^ Im(Aj (g bc + g)) Rc(A 3 (J bc + B)) 
^ \\ 3 {D hc + B)\± 



= 0, (1.18) 
etc. 

In Section 2 we succinctly consider abstract damped wave equations and de- 
tail the intimate spectral connections between (abstract versions of) GT bc . a / P 2 and 
D bc + B, a topic discussed in depth in [39]. Section 3 is devoted to the self-adjoint 
supersymmetric Dirac-type operator Z? bc + B in the Hilbert space L 2 ([0, 1]; p 2 dx) 2 
and the associated operators T hc . In particular, we describe in detail the boundary 
conditions chosen at x = 0, 1 and the Green's function corresponding to T£ T b , 
that is, the integral kernel of the resolvent (T bc T bc — zl)^ 1 . In Section 4 we derive 
an explicit formula for (an abstract version of) the resolvent of Dbc + B, employing 
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the well-known expression of the resolvent of the supcrsymmctric Dirc-type opera- 
tor Dbc- Section 5 focuses on the infinite sequence of trace formulas for the damped 
string equation (cf. (1.13)— (1.18)), using several well-known results in the spectral 
theory of non-self-adjoint operators associated with the names of Schur, Livsic, and 
Keldysh. In addition, based on abstract results on the existence of a Riesz basis 
with parentheses for a certain class of non-self-adjoint perturbations of a normal op- 
erator due to Katsnelson [62], [63], Markus [67], and Markus and Matsaev [69], [70], 
we show that Z?bc + B possesses a Riesz basis with parentheses in L 2 ([0, 1]; p 2 dx) 2 
(cf. Theorem 5.14). A host of useful (spectral) properties of abstract self-adjoint 
supersymmetric Dirac-type operators is collected in Appendix A. 

Finally, we briefly summarize some of the notation used in this paper: Let % 
be a separable complex Hilbert space, (•, -)u the scalar product in % (linear in the 
second factor), and I-h the identity operator in H. Next, let T be a linear operator 
mapping (a subspace of) a Banach space into another, with dom(T), ran(T), and 
ker(T) denoting the domain, range, and kernel (i.e., null space) of T. The closure 
of a closable operator S is denoted by S. The spectrum, essential spectrum, point 
spectrum, discrete spectrum, and resolvent set of a closed linear operator in % will 
be denoted by cr(-), c ess (-), <r p (-), <7d(-), and p(-), respectively. The Banach spaces 
of bounded and compact linear operators in H are denoted by B(H) and Boo(H), 
respectively. Similarly, the Schatten-von Neumann (trace) ideals will subsequently 
be denoted by B p i^K), p G (0, oo). Analogous notation £(%!, % 2 ), £>oo(%i, etc., 
will be used for bounded, compact, etc., operators between two Hilbert spaces T-L\ 
and 7^2 • In addition, tr^(T) denotes the trace of a trace class operator T 6 B\{H) 
and det-n y k(In + S) represents the (modified) Fredholm determinant associated 
with an operator S G BkiW.), k e N (for k — 1 we omit the subscript 1 in det^(-)). 
Finally, Pa/1 denotes the orthogonal projection onto a closed, linear subspace Ai of 
"H, + denotes the direct (not necessarily orthogonal) sum in *H, and abbreviates 
the direct orthogonal sum in H. 

2. Abstract Damped Wave Equations 
We start with some abstract considerations modeling damped wave equations. 

Hypothesis 2.1. Let % be a complex separable Hilbert space. 

(i) Assume that A is a densely defined, closed operator in % satisfying 

A* A > el n (2.1) 

for some e > 0. 

(ii) Suppose that R e B(H). 

We note that various extensions of the condition R 6 B(TL) are possible, see for 
instance, [8], [9], [33, Sect. VI.3]. 

Since kcr(A*A) = kci(A) if A is densely defined and closed in H, assumption 
(2.1) implies that 

ker(A) = kcr(A*A) = {0}. (2.2) 

In addition, we recall that (2.1) implies that ran(A) is a closed linear subspace of 
H. 

Since we are interested in an abstract version of the damped wave equation of 
the form 

u tt + Ru t + A* An = 0, u(0) = / , u t (0) = f u t > 0, (2.3) 
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we rewrite it in the familiar first-order form 



MM fu(0)\_ff 



K u t )-\-A*A -R)\u t )> \uM)~\h)' '- ' (2 ' 4) 

and set up the abstract initial value problem (2.3) as follows: First, one introduces 
the Hilbert space Ha by 

H A = {dom(A);(;-) nA ), (u, v) Ha = {Au, Av) H , u, v G dom(A) (2.5) 

(Ha is complete since by hypothesis, A* A > el-u)- Moreover, since ker(A) = {0}, 
the polar decomposition of A in H is of the form A = Ja\A\, where J a is an isometry 
in H and \A\ = {A* A) 1 ' 2 > e x l 2 I u . Thus, \\Af\\ n = \\\A\\\ n , f G dom(A), and 
hence one actually has 

H A = H lA \ with \A\ >s^ 2 I n . (2.6) 

We emphasize, that while we assumed that ker(A) = {0}, it may happen that 
ker(^4*) ^ {0} and hence ran(^4) C % ; since H = kcr(A*) © ran(A). In particular, 
we are not assuming G p(A), but (2.1) implies G p(|^4|). 

Given Ha, one then studies the abstract Cauchy problem in the Hilbert space 
Ha © H, 

F'(t) = G AiR F(t), F(t) G dom(G AiiJ ), t > 0, F(0) = (fy , (2.7) 

where 

Gam = ( Jj* M lH ^j , dom(G A , R ) = dom(A*A) © dom(A) CH A ®H, (2.8) 
= (^j) e dom(GU,*), t > 0, (2.9) 



and the scalar product in Ha © H is of course defined as usual by 

= {ui,u 2 )n A + ( v i, v 2)u 



Vl J ' I V2 



H A ®n 



= {Au 1 ,Au 2 ) n + {v u v 2 ) H , f^J ,(^ 2 2 J &H A ®H. (2.10) 

Since iGa,o is well-known to be self-adjoint in Ha®H, generating the corresponding 
unitary group in Ha © H, 

GAot ( cos ( y (A*A) 1 / 2 t) (A* A)- 1 / 2 sin ((A* Ay/ 2 t)\ 

~ \- (A* A) 1 / 2 sin ((A* Ay/ 2 t) cos ({A*A) 1 / 2 t) )' *~° 

(2.11) 

(cf. e.g., [33, Sect. VI.3], [44, Sect. 2.7], [45]), and 

B R = (° q EB(H A (BH), (2.12) 

Gam = Ga,o + Br generates a Co group on Ha © H (see, e.g., [33, Sect. VI.3], [44, 
Sect. 2.7], [45]). 

Next, we make a connection between the point spectral properties of iGam in 
Ha © H and an abstract perturbed Dirac-type operator Q + S in H © H defined by 

Q=\ A A q), dom(Q) = dom(A)ffidom( J 4*) C H®H, (2.13) 
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-iR 




dom(S) = H®H, (2.14) 



Q + S=(^ l A X)' dom (Q + S ) -dom(^)0dom(A*) C H © (2.15) 

For more details on the self-adjoint operator Q we refer to Appendix A. 
We start with the following elementary result: 



Lemma 2.2. Assume Hypothesis 2.1. Then 

-iR A*\ ( I n \ _ ( M./ /„ 
i4 y ^-zA J \-iA y ^-AM -R 

on dom(^*A) © dom(A) C "Ha © "H. 

Proof. A direct calculation yields for u G dom(A*A) and v G dom(^4), 

-£R A*\ ( 7«\ AA _ f-iA*Au - iRv 
A j I -a J \vj ~ \ Av 



and 



In] . ( \ /u\ _ / -iA*Au - iRv 

-iA j * l-AM -i?J U J ~ l Av 



(2.16) 

(2.17) 

(2.18) 
□ 



In this context one observes that the facts A* A = \A\ 2 and dom(^4) = dom(|A|) 
permit one to replace the Dirac-type operator Q by 

Q\A\=(\\\ 'o'), dom(Q| A |) =dom(A)©dom(yl) QH®H, (2.19) 

which can be advantageous as \A\ > e 1 / 2 /^ (cf. (2.6)). 

Next, we recall a few facts regarding eigenvalues of a densely defined, closed, 
linear operator T in H: The geometric multiplicity, m g (\ ,T), of an eigenvalue 
A G o-p{T) °f T is given by 

m g (\ ,T) = dim(ker(T- A /«)), (2.20) 

with ker(T — XqI-h) a closed linear subspace in T-L. (Here, and in the remainder of 
this paper, dimension always refers to the cardinality of an orthonormal basis (i.e., 
a complete orthonormal sequence of elements) in the separable pre-Hilbert space 
in question, cf. [104, Sect. 3.3].) 

The set of all root vectors of T (i.e., eigenvectors and generalized eigenvectors, 
or associated eigenvectors) corresponding to A G cr p (T) is given by 

K(\ ,T) = {/ £H \ {T-\ I u ) k f = Ofor some k G N}. (2.21) 

The set 1Z(Ao,T) is a linear subspace of H whose dimension equals the algebraic 
multiplicity, m a (Xo,T), of Ao, 

m a (X ,T) = dim ({/ eH\(T- \ {) I H ) k f = for some k G N}). (2.22) 

In general, 1Z(\q,T) is not a closed linear subspace of %. (1Z(\q,T) is of course 
closed if m a (Xo,T) < oo.) One has 

m g (X ,T) < m a (X ,T). (2.23) 
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If in addition, the eigenvalue Ao of T is an isolated point in <x(T) (i.e., separated 
from the remainder of the spectrum of T), one can introduce the Riesz projection, 
P(Xo,T) of T corresponding to A , by 

P(A ,T) = --^/ dCCr-C/w)- 1 , (2.24) 

2™ Jc(\„;e) 

with C(Ao;e) a counterclockwise oriented circle centered at Ao with sufficiently 
small radius e > 0, such that the closed disk with center Ao and radius e excludes 
cr(T)\{Ao}. In this case, 

H = kcr(P(A ,T))+ran(P(A ,T)), (2.25) 

with ker(P(A ,T)) and ran(P(Ao,T)) closed linear subspaces in H, and (cf. [80, 
No. 149]) 

ran(P(A ,T)) = {/e^| lim ||(T - X Q In) n f\\ 1/n = 0}, (2.26) 

and hence, 

n(\ ,T) Cran(P(A ,T)). (2.27) 

Moreover, in the particular case where 

m a (Xo,T) < 00, or equivalently, dim(ran(P(Ao, T)) < 00, (2.28) 

one has (cf. [42, Sect. II. 1], [43, Sects. 1.1, 1.2]) 

™ a (A , T) = dim(ran(P(A , T)) < 00, 1Z(\ ,T) = ran(P(A , T)). (2.29) 

If T G Boo(H), then any A G er p (T)\{0} satishes (2.29). 

We also note for later purpose that if A = {Ai, . . . , Ajv} for some TV G N, repre- 
sents a finite cluster of eigenvalues of T, isolated from the remaining spectrum of 
T, then the Riesz projection of T corresponding to A is given by 

P(A,T) = ~<f d((T-(I H )-\ (2.30) 

2m JC(A) 

with C(A) a counterclockwise oriented simple contour enclosing all eigenvalues in 
the cluster A in its open interior, such that C(A) excludes a(T)\A. In this situation 
P(A,T) represents the sum of the disjoint Riesz projections corresponding to Aj G 
A../ I V. 

iV 

P(A, T) = P(^j,T), P(Xj,T)P(Xk,T) = 6 j!k P(X k ,T), j, k = 1, . . . , N. 

(2.31) 

The connection between the point spectra of IGa,r and (Q + 5), more precisely, 
between %Ga,r and (Q + S)(I-u © P ra n(A))i is detailed in the following result. In 
this context one observes that Hypothesis 2.1 implies that ran(A) is a closed linear 
subspace of %, 

ran(^) = ran(A) = {kci{A*))^ . (2.32) 
Theorem 2.3. Assume Hypothesis 2.1. T/ien 

i <T p (iG A ,R), (2.33) 

and 

T P («i ; ii) = ( 7 P (O + 5)\{0}, (2.34) 
wzi/i geometric and algebraic multiplicities preserved. 
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More precisely, let ^ Ao G o~ v {iGA,R) and 

-*(«£>)■ (SO- 1 **- ,2 - 35 » 

i/ien A G cr p (Q + S) and 

(q + 5) Ll^J - ao Li^J . f j&U * + s) . ^ 



-iAu(X Q )J °\-iAu(\o))' \-iAu(Xo) 
Conversely, let ^ Ai G cr p (Q + S) and 

(Q + S) f^S^ = Ai ft 1 !. 1 !") . (ifA ) cbn f (0 4 .S... (2.87) 



then Ai G o- p {iGa,r), 4>2(Xi) G ran(A), and 

^"^ran^^CAl)^ _ . (lA~ 1 | r an(yl)^2 (Al ) 



(2.38) 

G dom(G^ iK ), 



V>i(Ai) 

where ^4 | ra n(A) denotes the inverse of the map A : dom(^4) — > ran(A). 

We note that Q + S in (2.34), (2.36), and (2.37) can be replaced by (Q + S)(I-u 

-fran(A))- 

Proof. Let (u(A ) w(A )) T G dom(GA,R) and suppose that iGA,ii(u(\o) w(A )) T 
0, that is, 



I n \ fu(X )\ ( v(Xo) 
-A* A -R) W(A W \-A*Au(X ) - Rv(X ) 



0, (2.39) 



implying u(X ) = v(X ) = 0, and hence (2.33). 

Next, let ^ Ao G o- p (iGa,r) and suppose (2.35) holds. Then, (u(Ao) u(Ao)) T G 
dom(GA,fl) implies u(Ao) G dom(^4*^4) and v(X ) G dom(^4) and using (2.16) one 
obtains 

-iR A*\ ( v(X ) \ _ (-iR A*\ (0 l n \ fu(X o y 
A oJ\-iAu(X )J \ A 0j\A J \v(X ) 

-iR A*\J I n \fu(X )\ =x J v(X ) (24() . 



A J \-A*A -Rj\v(X )J ""\-iAu{X ); 

proving (2.36). It is also clear that if (uj(X ) Vj(X )) T , j = 1, 2, are two linearly in- 
dependent nonzero solutions of iG a,r{u{Xq) w(Ao)) T = Ao(u(Ao) u(Ao)) T , then also 
( Wj (A ) - iA Uj {X )) T , j = 1, 2, of (Q + S)(M^o) MM) T = Ao(^i(Ao) ^ 2 (A )) T 
are linearly independent nonzero solutions since by hypothesis, ker(A) = {0}, prov- 
ing that geometric multiplicities are preserved. 

Finally, supose that ^ A x G cr p (Q + S) and assume that (2.37) holds. Then 
one obtains 

- tityi(Ai) + A*V 2 (Ai) = AiVi(Ai), ^i(Ai) G dom(A), V> 2 (Ai) G dom(A*), 

(2.41) 

A0i(Ai) = Ai^ 2 (Ai), implying ^(Ai) G ran(A) = dom^ 1 | ran(A) ) . (2.42) 
Thus, one computes 

o i H \ (iA-\ &Tl[A) M^T 

-A* A -R I Vi(Ai) 
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lA- 1 ^^ f-iR A*\ /Vi(Ai 

I H J\ A J y 2 (K 

^ _1 |ran(A)\ flpl(M)\ _ x ( iA~ 1 | ran (A) 1p2 (A 



Again one infers that geometric multiplicities are preserved as the map ^4 _1 | r an(A) : 
ran(A) — > dom(A) is injective. 

Finally, the preservation of algebraic multiplicities follows from the unitary equiv- 
alence result in Theorem 2.4 below. □ 

In fact, one can prove the following extension of Theorem 2.3, and we refer to 
[39] for a detailed proof: 

Theorem 2.4 ([39]). Assume Hypothesis 2.1. Then % © (kcr(A*))- L = "Hffiran(A) 
is a reducing subspace for Q + S and 



(Q + S)(I n © [I n - P ker{ A*)}) = (Q + S)(I n © PrMA)) 
= U x iG A , R U~\ 



(2.44) 



where A defined by 



t \Ua^ rant A), . . . . 

A: { v h is unitary, 2.45 

/ ^ Af, 



and 



U A=[_° i X 'o ) W ' W i-iin(.l).) w ».;//7 ( /; 7/ . (2. 46) 



U~ 1 = ^° M ^ G B(W © ran(A), H A © W) «s unitary. (2.47) 

We note that while Theorem 2.4 has been obtained by Huang [57] under the 
additional assumption \ A\ > e 1 / 2 /^, the results in Theorems 2.3 and 2.4 appear to 
be new under the general Hypothesis 2.1. 

Next, let £ be a conjugation operator in H, that is, £ is an antilinear involution 
satisfying (see, e.g., [30, Sect. III. 5] and [41, p. 76]) 

{£u,v)n = (£v,u)n u,v eH, € 2 = I n . (2.48) 

In particular, 

{£u,£v)n = (v,u)n, u,veH. (2.49) 
The densely defined operator S in H is called ^-invariant if 

S=<LS£. (2.50) 

Lemma 2.5. Assume Hypothesis 2.1, let £ be a conjugation operator in H and 
suppose that A* A and R are ^-invariant. Then 

\q G o- p (iG A ,R) if and only if - A G o- p (iGA,n), ( 2 -51) 

with geometric and algebraic multiplicities preserved. As a consequence of Theorem 
2A, the analogous relation (2.51) (including preservation of geometric and algebraic 
multiplicities) extends to nonzero eigenvalues of Q + S. 
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Proof. This follows from 

€)^^ H \0 £ 



" N '( '.\.n ") = /r;.,.„. (2.52) 

□ 



The following remark illustrates the connection between the generator iGa,r and 
the monic second-degree operator polynomial L(z) = z 2 I-u + ziR — A* A, z 6 C, in 
H: 

Remark 2.6. The eigenvalue problem for iGa,r, that is, 



U(Ao)J = Ao [v(Zy ' [v(Zy e doMGA ^ (2 - 53) 

is equivalent to the pair of equations 

v(X ) = -i\ u(X ), A*Au(X ) + Rv(X ) - iX v(X ) = 0. (2.54) 

Thus, eliminating the component v(Xq) in the second equation in (2.54) yields the 
quadratic pencil equation for w(A ), 

(A* A - X a iR - Xl)u(X a ) = 0, u(X ) f= dom(AM). (2.55) 

Similarly, eliminating ip2(Xi) in the eigenvalue equation (2.37) for (Q + S) results 
in the analogous quadratic pencil equation for ipi(X\), 

(A* A - XiiR - Ai)V'i(Ai) = 0, ^i(Ai) e dom(A*A). (2.56) 

Moreover, introducing the quadratic pencil (i.e., the second-degree operator poly- 
nomial) L(-) in %, 

L{z)=z 2 I H + ziR-A*A, dom(L(z)) = dom{A*A) 1 zeC, (2.57) 
one verifies the identity 

(L(z) © In)F(z) = E(z)(iG A ,R - zI H(Bn ), z e C, (2.58) 
where (for z e C) 

E{z) = (^ ZIh Ih lR ~ Q 7 ^) , dom(E(*)) - dom(.R) ®HCH®H, (2.59) 

fiW" 1 ^^ dom^z)- 1 ) = H®dom(R)CH®H, 

(2.60) 

Thus, identity (2.58) exhibits iGa,r as a global linearization of the quadratic pencil 
£(•) (in analogy to the discussion in [81, Sect. 1.1, Example 1.1.4] in the context of 
bounded operator pencils) , and again the (geometric and algebraic) multiplicity of 
nonzero eigenvalues of L(-) and %G a,r coincide by definition. Of course, the unitary 
equivalence described in Theorem 2.4 also exhibits (Q + S) ( I-u © -P ra n(A) ) as a global 
linearization of L(-). 

We note that even though the pencil L(-) has unbounded coefficients, replacing 
L(-) by L(-)(A*A + In) 1 reduces matters to a pencil with bounded coefficients, 
in particular, [68, Lemmas 20.1 and 20.2] apply to the spectrum of L(-) in this 
context. 
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Finally, we note that the standard separation of variables argument, making the 
ansatz 

u(t) - e- lXt u{\), t G R, A G C, (2.62) 
and inserting it into the equation 

;«),=M:'S)- ^ 

then yields of course the familiar eigenvalue problem 

-*(-$*>)' (2 ' 64) 

compatible with (2.53) and (2.54). 

3. Supers ymmetric Dirac-Type Operators 

In this section we study self-adjoint supcrsymmetric Dirac-type operators D of 
the type 

D= (t ^0*) ' dom ( £) ) = dom ( r ) © dom(T*) (3.1) 
in the Hilbcrt space 

L 2 ([0, 1]; p 2 dx) 2 = L 2 ([0, 1]; p 2 dx) ® L 2 ([0, l];p 2 dx) ~ L 2 ([0, 1]; p 2 dz) <8> C 2 , (3.2) 

where T (and its adjoint T*) are densely defined and closed first-order differential 
operators in £ 2 ([0, l};p 2 dx) of the form 

T=-± T* = \±p, (3.3) 
p dx p l dx 

with appropriate boundary conditions at x = 0,1, and p > is an appropriate 
weight function, which, throughout this paper, will be assumed to satisfy the fol- 
lowing conditions: 

Hypothesis 3.1. Suppose that 

< p G Z/°°([0, 1]; dx), 1/p G L°°([0, 1]; dx). (3.4) 

Assuming Hypothesis 3.1, we now introduce the following concrete models for 
the operator T in L 2 ([0, 1]; p 2 dx): 

T max / = (i/p)f, (3.5) 

/ G dom(T max ) = {g G i 2 ([0, l];p 2 dx) \ g G AC {[0,1]); g' G L 2 ([0, l];p 2 dx)}, 

T mia f = (i/p)f, 

f G dom(T min ) = {g G L 2 ([0, l];p 2 dx) \ g G AC([0, 1]); g(0) = = g(l); (3.6) 

g' eL 2 ({0,1}; p 2 dx)}, 

Tof = (i/p)/', 

/ G dom(To) = {g G L 2 ([0, 1]; p 2 dx) \ g G AC([0, 1]); 3 (0) = 0; (3.7) 

3 'Gi 2 ([0,l];p 2 dx)}, 

Ti/ - (i/p)f, 

f G dom(T 1 ) = {g G L 2 ([0, 1]; p 2 dx) \ g G AC([0, 1]); g(l) = 0; (3.8) 

g' &L 2 ({0,1}; p 2 dx)}, 
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T u f = (i/p)f, ^C\{0}, 

/ G dom(T w ) = {g e L 2 ([Q, l];p 2 dx) \ g G AC([0, 1]); g(l) = ug(0); (3.9) 

g' EL 2 ([0,1]; p 2 dx)}, 

where AC([0, 1]) denotes the set of absolutely continuous functions on [0, 1]. Due 
to the conditions (3.4) on p, the fact that all operators in (3.5)— (3.8) are closed and 
densely defined in L 2 ([0, l};p 2 dx) parallels the well-known special case where p = 1 
a.e. 

The associated adjoint operators to (3.5)-(3.8) are then given by 
T^J = (i/p 2 )(pf)', 

f G dom(T* ax ) - {g G L 2 ([0, l];p 2 dx) \ pg G AC([0, 1]); (pg)(0) = = (pg)(l); 

(pg)' EL 2 ([0,1]; p 2 dx)}, (3.10) 

Tmin/ = (VP 2 )(P/)', (3.11) 

/ £ dom(T* in ) = {g G L 2 ([0, 1]; p 2 dx) \ pg E AC ([0,1]); (pg)' G L 2 ([0, 1]; p 2 dx)}, 
T ~f = (i/p 2 )(pf)', 

f G dom(T*) = {g G L 2 ([0, 1];^) | Pff G AC([0, 1]); (p<?)(l) = 0; (3.12) 

(pg)' EL 2 ([0,l];p 2 dx)}, 

nf = (i/p 2 )(pfy, 

f G dom(TT) = {g G L 2 ([0, l];p 2 cte) | A9 G AC([0, 1]); (p<?)(0) - 0; (3.13) 

(pg)' eL 2 ([0,l];p 2 dx)}, 

T*J = (*/p 2 )(pf)', w G C\{0}, 

/ G domCO = {.9 G L 2 ([0, l];p 2 dx) \ pg G AC([0, 1]); (pg)(l) = (l/uj)(pg)(0); 

(pg)' eL 2 ([0,l};p 2 dx)}. (3.14) 

The sufficiency of the boundary conditions in (3.10)-(3.14) is clear from an ele- 
mentary integration by parts 

[ p(x) 2 dxj(x)[(l/p(x))g'(x)] = f dx[p(x)f(x)]g'(x) 
Jo Jo 

t-1 



= [p(x)f(x)]g(x)y - / dx[p(x)f(x))'g(x) 
Jo 

= [p(x)f(x)]g(x)\l- [ p(x) 2 dx[l/p(x) 2 ][p(x)f(x)]'g(x). (3.15) 
Jo 

We omit the details of the necessity of these boundary conditions which are well- 
known in the special case p = 1 a.e. 

Replacing T in (3.1) by T max , T m i n , T , T\, and T u , the corresponding self-adjoint 
Dirac-type operators in £ 2 ([0, 1]; p 2 dx) 2 are then denoted by D max , -Dmin, Do, D\, 
and D w , respectively. 

Next we analyze the nullspaccs of all these operators: 

Lemma 3.2. Assume Hypothesis 3.1. Then 

kcr(T max ) = C, kcr(T* ax ) = {0}, (3.16) 
ker(T min ) = {0}, kcr(T* in ) = {c/p | c G C}, (3.17) 



THE DAMPED STRING PROBLEM 



13 



and hence 



kcr(T o ) = {0}, ker(T *) = {0}, (3.18) 
ker(Ti) = {0}, ker(T 1 *) - {0}, (3.19) 

ker(Tw) = ({0}, ^C\{0,1}, = f{0}, ^C\{0,1}, 

\C, w = l, " \{c/ P \ceC}, w = l, 

(3.20) 

ker(.D max ) = C © {0}, (3.21) 
ker(£ min ) = {0} © {c/p | c G C}, (3.22) 
ker(A)) = {0}, (3.23) 
ker(Di) = {0}, (3.24) 

ker(D w ) = (j?*' , ^ ^C\{0,1}, 

IC©{c/p|ce C}, w = l. 



Proof. Since 



d cl=0. 



p(x) \dx J p(x) 



d 



0, ceC, (3.26) 



dx\ p{X \p(x), 

in the sense of distributions on (0, 1), and in each of these cases the operators 
T and T* are of first order and hence have at most a one-dimensional nullspace, 
one concludes that C C L 2 ([0, 1]; p 2 dx) (resp., {c/p \ c G C} C £ 2 ([0, 1]; p 2 ^)) 
are the kernels for T (rcsp., T*) if and only if c (resp., c/p) satisfy the boundary 
conditions in T (resp., T*). Checking whether the boundary conditions are fulfilled 
is elementary and yields (3.16)-(3.20). 

Relations (3.21)-(3.25) then follow from the general fact that ker(D) = ker(T) © 
ker(T*) (cf. (A.23)). □ 

We continue by listing the corresponding operators T*T: 

r* ax r max / = -(l/Af, 

/ G dom(T* ax T max ) = {g G L 2 ([0, l];p 2 dx) \g,g'e AC([0, 1]); (3.27) 
g'(0) = = g l (l);g" &L 2 ([0,l];p 2 dx)}, 

r* in T min / = -(i/p 2 )/", 

/ G dom(T* in T min ) = {g G L 2 ([0, l];p 2 dx) \g,g'e AC([0, 1]); (3.28) 
g(0) = = g(l);g" GL 2 ([0,1]; p 2 dx)}, 

T*T f = -(l/p 2 )f", 

f G dom(T*T ) = {g G L 2 ([0, 1]; p 2 dx) \ g, g' G AC ([0,1]); (3.29) 

g(0) = = g'(l);g" e L 2 ([0,1]; p 2 dx)} , 
T*TJ = -(l/p 2 )f, 

f G dom^TJ = {. 9 G L 2 ([0,l];p 2 dx) \g,g' G AC([0,1]); (3.30) 
3 '(0)=0 = . 9 (l);.g"GL 2 ([0,l];p 2 dx)}, 



1 For simplicity, denotes the zero vector in L 2 ([0, 1]; p 2 dx) as well as in L 2 ([0, 1]; p 2 dx) 2 . 
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t:tj = -(i/p 2 )/", 

/ G domCCTJ = {geL 2 ([0,l];p 2 dx)\g,g' GAC([0,1]); (3.31) 
g(l) = wg(0), g'(l) = (l/^'(0); 5 " G L 2 ([0, 1]; p 2 dx)}. 
Analogously, one obtains for the operators TT*: 

T m ^j = -(i/p)[(i/p 2 )( P fy]', 

f e dom(T* ax T max ) = { 3 G L 2 ([0, l];p 2 <fe) | pg,p- 2 (pg)' G AC([0, 1]); (3.32) 
(pg)(0) = = (pg)(l); [(1/ p 2 )(pg)'}' G L 2 ([0, 1]; p 2 dx)} 1 

T mia T^j = -(i/p)[(i/p 2 )(pfy}', 

f G dom(T min T min ) - {g G L 2 ([0, 1]; p 2 dx) \ pg,p- 2 {pg)' G 4C([0, 1]); (3.33) 

(pg)'(0) = = (pg)'(l); [{I / p 2 ){pg)']' G i 2 ([0, l];p 2 dx)}, 
T T*f = -(l/p)[(l/p 2 )(pfY}', 

f G dom(T*T ) = {g G L 2 ([0, 1]; p 2 dx) | pg, p- 2 (pg)' G AC([0, 1]); (3.34) 

(pg)'(0) = = (pg)(l); [{I / p 2 ){pg)']' G L 2 ([0, l]; P 2 dx)}, 
T 1 T*f = -(l/p)[(l/p 2 )(pf)'}', 

f G dom^Ti) = {g G L 2 ([0, 1]; p 2 dx) \ pg, p- 2 {pg)' G AC([0, 1]); (3.35) 
(p</)(0) = = (pg)'(l); [(l/p 2 )(pg)'}' G L 2 ([0, l];p 2 cfe)}, 

t w t:/ = -(i/ p )[(i/p 2 )(p/)']', 

/ G domCCTJ = {.9 G L 2 ([0, 1]; p 2 dx) \ pg, p- 2 (pg)' G AC([0, 1]); (3.36) 

( W )(l) = (l/5J)(w)(0), (pff)'(l) = w( W )'(l); [(l/p 2 )( Pff )Tei 2 ([0,l];p 2 dx)}. 

One recalls that the spectra of all these operators T*T and TT* in (3.27)-(3.36) 
are purely discrete, that is, 

^(^max,min,0, l,w^max,min,0,l,Lj ) ^"d (-^~rnax,min,0, l,w^max,min,0, l,w ) ' 

CT ess(^max,min,0,l,w^max,min,0,l,w) = $J (3.37) 

^(-^max,min,0, l,w~^max,min,0,l,w ) = M T max,minAi,w T max,min,o,i,u>)' 

CT ess(^max,min,0,l,w^max,min,0,l,w) = 0; (3.38) 

in addition, we have of course (cf. (A. 8) and (A. 10)) 

^(^max,min,0.1,Lj-^max,min,0,l,a;)\{0} ^(^max, min,O,l,w^max,min,O,l,w)\{0} (3.39) 

and 

kcr(T*T) = ker(T), kcr(TT*) = kcr(T*). (3.40) 
Here, and occasionally later on, T stands for T max , T m j n , T , T\, or T u . 

For subsequent purpose we next describe the connection between T m ax,min,o,i,w 
and T max min 01(J in the Hilbert space £ 2 ([0, 1]; p 2 dx) and the corresponding oper- 
ators 

i-j-) = r ma x,min,o,i,w(p = 1) in L 2 ([0, l];dx), 

"' X / max, min, 0,1, oi . . 

(d \ * ' 
* ) = T max,min,0,l,u,(/> = *) in ^(I , ^ ^)) 
a2; / max,min,0,l,w 
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taking p = 1 in (3.5)-(3.14). Introducing the unitary operator U p via 
. iL 2 ([0,l];p 2 dx)^L 2 ([0,l};dx) 



one obtains 



\ dX J max,min,0,l,w 



(3.43) 



^ P max, mm 



max,min,0.1.Lj 



where M w denotes the bounded operator of multiplication by the function w <E 
L°°([0, l];dx) in the Hilbert space £ 2 ([0, 1]; dx), and hence, 

/ (f \ 

^^max,min,0,l,w^max,min,0,l,w^CJ = -^1/p I ~ ~T~a ) ^1/p; (3-44) 

V "2- / max,min,0,l,w 



^P^max,min,0,l,w^max,min, 0,l,u)^f> ~ * j ) ^/P 2 ( * J™ 



dx I \ dx 

' ' max,min,0,l,o; V / max,min,0,l,u> 



where we abbreviated 

d 2 \ (A V ( . d 



dx 2 ) V dx / \ dx 

/ max, min, 0,1, w V / max,min,0,l,w \ / max, mill, 0,1, w 

= ^max,min,0,l,w(P = 1 )2niax,min,0,l,w (P = 1) (3.45) 

in L 2 ([0, 1]; dx), which are obtained as in (3.27)-(3.31), taking p = 1. 
Both relations in (3.43) immediately follow from a consideration of 

(f,Tg) L 2 mi] . p 2 dx) , f EL 2 ({0,1}; P 2 dx), gedom(T) (3.46) 

respectively, 

(/,T* 5 ) i2([0>1] . p2(te) , /6l 2 ([0,l];A), sedom(r). (3.47) 
Relation (3.44) implies the following trace ideal result: 
Theorem 3.3. Assume Hypothesis 3.1 and let 

^ ^- /^(-^max, mill, 0,l,w^niax, mill, 0,1, w ,min,0, l,W"^max,min,0, 1,oj 

). (3.48) 

Then 

(^max,min,0,l,w^max,min,0,l,a; — ^ ^l(-^ 2 ([0, l];p 2 rfx)), 

(^max,min,0,l,w^max,min,0,l,w — £ ^l(^ 2 ([0, l]',p 2 dx)). 

and 

7max,min,0,l,o;(T maXim i nj o,l,w^max,miii,0,l,w — Z "0 ^ ^2 ( [0, 1] J p 2 dx)) , 
^max,min,0,l,a;(^max,min,0,l,w^max,min,0,l,w — Z "0 ^ ^2 ( [0, 1] 5 p 2 dx)) . 

Proof. Since (3.49), in the case of T^ ax min 01 w T max min l w , is well-known if p = 1 
(in which case one can explicitly determine the asymptotics of the nth eigenvalue 
involved as 0(n 2 ) as n — > oo), (3.49) follows from (3.44) and the fact that the 
nonzero eigenvalues of (TT) 1 / 2 and (TT*) 1 / 2 (i.e., the singular values of T and T*, 



(3.49) 



(3.50) 
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respectively) coincide, as recalled in (3.39). Using the familiar polar decomposition 
formulas 

T = V T \T\, T* = (V T )*\T*\, where \T\ = (T*T) 1/2 , \T*\ = (TT*) 1/2 , (3.51) 

for any closed densely defined operator T in H, with Vr a partial isometry in H 
(cf. [61, Sect. VI.2.7]), one obtains 

T(T*T - zIuV 1 = V T \T\(\T\ 2 - zI H )~\ 

1 (o.bZ) 

T*(TT* - zlu)- 1 = (V T )*\T*\(\T*\ 2 - zl n ) , 
proving (3.50). □ 

Of course, (3.49) and (3.50) extend to all z in the corresponding resolvent set of 
T*T and TT* , respectively. 

We conclude this section by listing the Green's functions (i.e., the integral kernels 
of the resolvent of) T*T at z — whenever the corresponding nullspaces are trivial 
(cf. (3.16)-(3.20)). We write 2 

((T*T - zl)-\f) (x) = f 1 p(x') 2 dx' G T , T (z, x, x')f(x'), (3.53) 
Jo 

ze P (T*T), f€L 2 ([0,l];p 2 dx), 

where T represents one of T m i ni0 ,i, w , and record the following result: 

Lemma 3.4. Assume Hypothesis 3.1. Then, 

• , \x(l-x'), 0<x<x'<l, . . 

G T *. t . (0,x,x') = { ' " - " ' 3.54 

G t *t (0, x,x') = \ X ] °^ x ^ x '^^ (3.55) 

Gt* Ti (P,x,x)= \ (3.56) 
\{l—x), < x < x < 1, 



, 1 , 1 + w-Lj- \uj\ 2 I-cj + cj- H 2 , 
Gt , t J0,x,x) = --\x-x\ x x 



Proof. In all cases depicted in (3.54)-(3.57), one has (T*T)- X e B(L 2 ([0, l];p 2 dx)) 
and hence Gt*t(0,x,x') is well-defined. Since T*T is of the type -(1/ p 2 )(d 2 /dx 2 ), 
one infers that for all c, d <E C, T*T(c + dx) = is valid in the sense of distributions, 
and hence (3.54)-(3.56) readily follow from the separated boundary conditions im- 
posed in (3.28), (3.29), and (3.30). The nonseparated boundary conditions in T*T U 
in (3.31) require a slightly different strategy: Introducing 



g u (x;f)= J p(x') 2 dx'GTzTjO,x,x')f(x'), f e L 2 ([0, 1]; p 2 dx), (3. 



58) 



We denote by I the identity operator in L ([0, 1]; p dx)) and similarly, by 1% the identity 
operator in L 2 ([0, 1]; p 2 dx) 2 . 
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g(x;f) = ~^p(x') 2 dx'\x-x'\f(x'), f G L 2 ([0, 1]; p 2 dx), (3.59) 
one infers that 

-pW^ [gu{x ' J) ~ 9{x ' J)] = 0, (3 - 60) 

and hence, 

g u (x;f)=g(x;f) + c 1 (f)+c 2 (f)x, ie[0,l], (3.61) 

for some coefficients Cfc(/) G C, fc = 1, 2. 

Since by definition, g u (-; /) G dom(T*T w ), the boundary conditions (cf. (3.31)) 

g u {l;f)=wg u {0;f), <£(1; /) = {l/ZJ)g' u {0; /), (3.62) 
and Cramer's rule determine Ck(.f), k = 1,2, implying (3.57). □ 

4. Non-Self- Adjoint Perturbations of Supersymmetric 
Dirac-Type Operators 

In this section we present an abstract result 3 , the computation of the resolvent 
of a non-self-adjoint perturbation of special diagonal type of a supersymmetric 
Dirac-type operator. 

Hypothesis 4.1. Le£ 'Hj, j = 1,2, be separable complex Hilbert spaces and let 

A : dom(A) C Hi -> H 2 (4.1) 

&e a densely defined closed linear operator. In addition, assume that 

V G (4.2) 

Given Hypothesis 4.1, we define the self-adjoint supersymmetric Dirac-type op- 
erator Q in Hi ® %2 by 

= ( °a A n\ dom(Q) = dom(A) ® dom(A*), (4.3) 



< A , 

introduce the special diagonal operator W G B(Hi © %) via 





and finally consider the perturbed Dirac-type operator 
Q + W = 



W=[ V n !), (4-4) 



G 










J - 







(4.5) 

dom(Q + W) = dom(Q) = dom(A) ® dom(A*). 
One then computes the following expression for the resolvent of Q + W: 
Theorem 4.2. Assume Hypothesis 4.1 and choose £ G C swc/i that 

C 2 E p(A*A)n P (AA*) and [l + (V(A*A - C^Mi) -1 ] ~* G B(H\). (4.6) 
Then Q G p(Q + and 

(g + ^-C/^e^)" 1 



3 The notation used in this abstract part is also employed in Appendix A. 



18 F. GESZTESY AND H. HOLDEN 

c(a*a-c 2 / Ki )- 1 -aA* A-ein 1 r 1 [i+c,v(A' A-ein 1 r 1 r 1 \ 

xlJ+CV^'yl-C 2 /^!)- 1 ]- 1 xVA'iAA'-^^-'+A'iAA'-ei^r 1 

AiA'A-^In,)- 1 -A(A* A-eiuJ-^I+CViA* A-ei^)- 1 }- 1 

. x[/+Cy(A*A-C 2 /w 1 )- 1 ]- 1 xVA*(AA*-C 2 I W2 )- 1 +C(AA*-C 2 I W2 )- 1 / 

Proof. A direct computation, using (A. 27) reveals 

/C^M-C 2 /^)" 1 A^AA*-^)" 1 

- \ y A(A*A - C 2 /^)" 1 C(AA* - C 2 /^)" 1 

7 Wl \ /V 0\ faA*A-C 2 I ni )- 1 A*(AA* -Clu,)- 1 
J W J + {0 Oj {A(A*A - C'lnJ- 1 C(AA* - C 2 ^)" 1 

fc(A*A-ei Hl )- 1 a*{aa* -eiu.r 1 
\A{A*A-ei Hl )- 1 c(aa* -ein.r 1 

I Hl \ (QV{A*A - ei^r 1 VA*{AA* - C 2 /^)" 1 

o /«J + V 

/C(A*A-C 2 /« 1 )- 1 A*(AA*-C 2 /« 2 )- r 

- {a(A*A - ein,)- 1 C(AA* - C 2 /^)- 1 

x / [/+CV(A*A-C 2 /« 1 )- 1 ]- 1 -[I+CViA'A-eiu^-^VA'iAA'-eiu^- 1 

V o / 



(4.7) 



(4.8) 
(4.9) 



using 

,0 InJ ~ { J Wa 
for F.F- 1 e B(Hi), G E B(H 2 ,Hi). Relation (4.8) proves (4.7). □ 

Of course this result extends to more general perturbations V: For instance, 
comparing \Q\ = ^ |^»| ^ and W, it is clear that V being relatively bounded with 

respect to (A* A) 1 / 2 with relative bound strictly less than 1, would be sufficient. In 
particular, the term 



VA*(AA* - C 2 ^)" 1 = [V(A*A - ( 2 I ni r 1/2 ] [(A* A - ( 2 I Hl )-V 2 A*] e B{Hi) 

(4.10) 

in (4.7) is then well-defined. More generally, one could invoke form rather than 
operator perturbations W of \Q\. 

5. An Infinite Sequence of Trace Formulas for 
the Damped String Equation 

The principal aim of this section is the derivation of an infinite sequence of trace 
formulas for the damped string equation. 

Throughout the major part of this section the coefficients a and p in (1.1) will 
be assumed to satisfy the following conditions: 

Hypothesis 5.1. Suppose that 

a e L°°([0, 1]; dx), a real-valued a.e. on (0,1), (5-1) 

and 

< p G L°°([0, 1]; dx), l/pe L°°([0, !];<&;). (5.2) 
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Introducing the sclf-adjoint Dirac-type operator D in L 2 ([0, l];p 2 dx) 2 by 

D= (t ^ ' dom (' D ) = dom ( T ) © dom(T*), (5.3) 

the concrete non-self-adjoint perturbations B of D then will be of the special diag- 
onal form 

B= (-ia/p 2 0j €B ( iJ([0)1 ] iW ). (5.4) 

Given Hypothesis 5.1, we introduce the closed (densely defined) Dirac-type op- 
erator 



T J V °J ' (5.5) 
dom(L> + B) = dom(£>) = dom(T) © dom(T*) 

in i 2 ([0, 1]; p 2 dx) 2 , where T represents one of T min , T , T\, and T u , as singled out 
in Lemma 3.4. 

Remark 5.2. To make the connection with the abstract damped wave equation 
discussed in Section 2, one identifies, % and L 2 ([0, l];p 2 dx) 2 , A and T, Q and D, 
B and S 1 , and R and —ia/p 2 , respectively. 

Our main aim is to derive an infinite sequence of trace formulas associated with 
D + B, but first we note the following result, denoting by 

hn(S) = (S-S*)/(2i), SeB(H) (5.6) 

the imaginary part of 5: 

Theorem 5.3. Assume Hypothesis 5.1 and choose (eC sucft i/iai 

C 2 6 p(T*T) n p(TT*) and [i + (B{T*T - C 2 /) -1 ] _1 £ B(Wi). (5.7) 
T/ien C G P(-D + B) and 
(D + B-Ch)- 1 



/ ((T*T-C 2 I)~ 1 II+(B{T*T-( 2 I)- 1 ]- 1 -C{T*T-( 2 I)- 1 [I+(B{T*T-C 2 I)~ 1 ]~ 1 
' y.BT*(TT*-C 2 iy 1 +T*(TT*-C 2 iy 1 

T(T"T-C 2 I)~ 1 [I+CB(T*T-C, 2 I)- 1 ]- 1 -T(T*T-C 2 /)~ 1 [/+CB(T*T-f 2 7)- 1 ]- 1 
\ xBT*(TT*-C 2 /) _1 + C(7 , T*-C 2 /) -1 



(5.8) 



with 

(D + B-Ch)- 1 £B 2 (L 2 ([0,1]; p 2 dx) 2 )\B 1 (L 2 ([0,1]; p 2 d:r) 2 ), (5.9) 

Imfp + B-C/2)- 1 ] eSi(L 2 ([0,l];p 2 da;) 2 ). (5.10) 

HereT and D stand for one o/T maX;m i n; o,i, w and Anax,min,o,i,w; respectively. 

Proof. Equation (5.8) and (D + B - Ch)' 1 S B 2 (L 2 ([0, 1]; ^dx) 2 ) are immedi- 
ate consequences of Theorems 3.3 and 4.2. To show that (D + B — Ch)^ 1 
Bi(L 2 ([0, l];p 2 dx) 2 ) we first focus on the term T(T*T - C 2 d) _1 in the (l,2)-entry 
of (5.8) and argue as follows: Since the eigenvalues of T*T are known to be asymp- 
totically of the form cn 2 + o(n 2 ) for some c ^ as n — > oo, one concludes that 

| T |( T * T + /)-i = |T|(|T| 2 + /) _1 i Si(L 2 ([0,l];p 2 da;) 2 ). (5.11) 

Next, we recall the polar decomposition for the densely defined and closed operator 
T, that is, T = V T \T\, or, |T| = (V T )*T, where \T\ = (TT) 1 / 2 > and V T is a 
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partial isomctry (cf. [61, Sect. VI.2.7]). Hence, if one argues by contradiction and 
assumes that T{T*T+I) _1 G £>i(i 2 ([0, 1]; p 2 dx) 2 ) , multiplication of \T\ (T*T+I) ~ 1 
from the left by the partial isometry (Vr)*, and applying the trace ideal property 
of Si(-), would yield the contradiction \T\(T*T + I)- 1 G #i(L 2 ([0, 1]; p 2 dx) 2 ). 
Replacing T by T* this also yields that T*(TT* + I)- 1 £ Si(L 2 ([0, l]-p 2 dx) 2 ) in 
the (2, l)-entry of (5.8). Since all other terms in (5.8) are in #i(L 2 ([0, 1]; p 2 dx) 2 ), 
this proves (5.9). 

Finally, (5.10) follows from (5.9) and the resolvent identity 

ImKD + B-Ch)- 1 ] =(D + B-(hy 1 [lm(C)-lm(B)}[(D + B-(h)- 1 Y. (5.12) 

□ 

Next, we turn to the expansion of tr^Qo iy p 2 dx )2 (im [(D + B — C-^) -1 ]) as 

Theorem 5.4. Assume Hypothesis 5.1 and choose ( G ^\{0} with \(\ sufficiently 
small, such that £ G p(D + B). Then, 

tr L 2([o,i];p 2 rfz) 2 ( Im [(D + B — (h)' 1 ] ) 

= Im [tr L2([0il] ., 2dx) ((2C + i(a/ p 2 ))(T*T - ( 2 I - (i(a/ p 2 ))" 1 )] 

OO 

= E^mC 2m , (5-13) 
m=0 

where 

to = ^^^^((a/^frr)- 1 ), (5.14) 

+ 3tr L2([0il] ., 2dx) ((a/p 2 ) (T*T)- 2 ) , (5.15) 
etc. 

HereT, D, andt m , m G N 0; stand for one o/T mini0 ,i, w , Anin,o,i,"; and t m in,o,i,w;m, 
to G No, respectively, and we only consider the case u> G C\{0, 1} (cf. (3.20)). 
Explicitly, one obtains for the coefficient t : 

£min ; o = / dx x(l - x)a(x) , (5.16) 

(5.17) 
(5.18) 

l )x + l]a(x), wgC\{0, 1}. (5.19) 

Proof. Using (5.8), one computes 

trL 2 ([o,i];P 2 ^) 2 ( Im [(D + B — Chy 1 ] ) 

= tv m%1] , p , dx) ( Im [C(T*T - C 2 /)- 1 [J - ^(a/p 2 ) (T*T fiy 1 ] ) 

+ tv LH[0Ahp2dx) (im [t(t*t - c 2 /)- 1 [J - C*Wp 2 )(r*T - c 2 !)- 1 ] - 1 



min;0 / 


(is x(l — x)a(x), 


Jo 




U):0 = 


dx xa(x), 


Jo 








*1j0 = / 


dx (1 — x)a(x), 


Jo 
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xi(a/p 2 )T*{TT* - fiy 1 
= tr L2 ([o,i];P^) (im [c(T*T - ( 2 I - Ci(a/p 2 )) _1 ] ) 
+ tr i2([0il];p2dx) ( Im [ T *T(T*T - C 2 /)" 1 

xiT-T-fl-CHa/f))- 1 ^//?) 
= tr^p,!];^) (im [C(T*T - C 2 / - C^//))" 1 ] ) 

- tr i2([0 , 1];p2dx) (lm [(T*T - C 2 / - C*(a/p 2 )) _1 H)(a/p 2 )]) 
+ tr i2([0il] .^ dx) (im [C 2 (T*T - C 2 /)" 1 

x{ T* T -C 2 I-^(a/p 2 ))-\( a lp 2 )\) 
= tr L 2([o,i];p 2 d») (im [(C + i(a/p 2 )) (T*T - C 2 / - C^/p 2 ))" 1 ] ) 
+ trx,2([ 0) i] ; p2 da .)^Im C^*? 1 - C 2 I - (i(a/p 2 )) 1 

xcila/p^rr-c 2 /)- 1 ]) 

= tr ia([ o,i ]ip2(fa) (im [(C + i(a/ p 2 ))(T*T - ( 2 I - Ci(a/p 2 )) _1 ] ) 

- tr L2([0 ,i ]; ,2 dx) (im [C(T*T - C 2 / - C^/p 2 ))" 1 
x (T*T - ( 2 I - (i(a/p 2 ) - T*T + ( 2 1)(T*T - ( 2 1)- 1 ) 

= tr ia([0 ,i] ;pa<fa ) (im [(C + i(a/ p 2 )){T*T - ( 2 I - ^(a/p 2 ))" 1 ] ) 

- tr i2([0jl];p2(ix) (im [C(T*T - C 2 /)" 1 - C(T*T - C 2 / - C^/p 2 ))" 1 ] ) 
tr L2([0 ,i];^) ((C + i( a / p 2 ))(T*T - ( 2 I - ^(a/p 2 ))- 1 )' 

}r L2mihp 2 dx) (C(T*T - C 2 /)- 1 - ((T*T - ( 2 I - (tia/p 2 ))- 1 )' 
tr i2([0 ,i];p^x) ((2C + i{a/p 2 )){T*T - ( 2 I - ^(a/p 2 ))" 1 ) 
= tr i2([0 ,i ];p2dx) (im [(2C + i(a/ p 2 ))(T*T - ( 2 I - ^(a/p 2 ))" 1 ] ) , (5.20) 

where we repeatedly used cyclicity of the trace (i.e., tr^ 2 (AB) = tr-n^BA) for A G 
B(n 1 ,H 2 ), B e B(H 2 ,Hi) with AB G Bi(W 2 ), #A G Bi(Wi), cf. [93, Corollary 
3.8]), (A.16), and the fact that tr w (Im(S)) = Im H (tr(S)) for S G Bi(ft). 

The following elementary computation, again employing cyclicity of the trace, 

Im 



= Im 

-Im 
= Im 



tr^ao.i];^*) ((2C + i(a/p 2 ))(T*T - ( 2 I - (i(a/p 2 )) ^ 
= tr i2([0 ,i ];p2<te ) (im [(2C + i(a/p 2 ))(T*T - ( 2 I - ^(a/p 2 ))" 1 ] ) 
- ^tr i2([0;1] ., 2dx) ((2C + t(a/p 2 ))(T*T - ( 2 I - &(<*/ p 2 ))" 1 

- (T*T - C 2 / + ^/p 2 )Y\2CI ^(a/p 2 ))) 
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- |tr i3 ([o,i];p"dx) ((T*T - C 2 / - C^/p 2 ))" 1 - - C 2 / + C^/p 2 ))" 1 ) 

+ \^m [Wdx) ({a/p 2 ) (T*T - C 2 / - (iia/p 2 ))- 1 

+ (T*T-C 2 I + ^(a/p 2 )y 1 (a/pi)) 
= j tr i2([0;1];p 2 dx) ((T*T - C 2 / - ^(a/p 2 ))- 1 - (T*T - ( 2 I + (i^/p 2 ))- 1 ) 
+ \tr LH[Wdx) ((a/p 2 ) (T*T - ( 2 I - (iia/p 2 ))- 1 

+ (a/p 2 ) (t*T ( 2 I + ^(a/p 2 ))- 1 ) , (5.21) 

then shows that the expression (5.20) is indeed even with respect to £ G R, proving 

(5.13) . 

Finally, expanding the leading terms in (5.13) up to order 0(( 2 ) then yields 
(5.15), and the expressions (5.16)-(5.19) are a consequence of (3.54)-(3.57) and 

(5.14) . □ 

Remark 5.5. Assume Hypothesis 5.1. 

(i) With respect to the location of the eigenvalues \j(D + B), j G J, we note 
the standard fact (cf. [61, Problem V.4.8]) that self-adjointness of D and B 6 
B(L 2 ([0, 1]; p 2 dx) 2 ) yield for d(z,a(D)) > ||5||L 2 ([o.i];p 2 d:r) 2 (with d(z, S) the dis- 
tance between z G C and the set £ c C) that 

\\B{D - zh)-^^^.^ < WBWv^.^idiz^iD))- 1 < 1, (5.22) 
and hence by the usual geometric series argument, 

\\( D + B ~ zI ?) 1 |lB(L2([o !l ]. p 2 da , ) 2 ) 

oo 

< d(z, a(D))~ 1 Y^iWBh^iy^Mz^iD))-^ 

k=0 

= [d(z, a(D)) - \\B\\ L 2 ([0tl] . p2dx) 2] ~\ (5.23) 
Thus, the (necessarily discrete) spectrum of D + B is contained in the strip, 

e(D + B) C {z G C | Im(z) G [ - \\B\\ L 2 ([QylWdx) 2, ||B|| L 2 ([0il]ip3<Ix) 2] }. (5.24) 

(ii) Moreover, since T*T and a/p 2 are invariant with respect to the operation of 
complex conjugation in L 2 ([0, l];p 2 dx) for T = T m i ni o,i, Lemma 2.5 applies to all 
nonzero eigenvalues of D + B and hence 

An e<7p{D + B) if and only if - A^ G cr p (L> + 5), A G C\{0}, (5.25) 

with geometric and algebraic multiplicities preserved. Here D represents one of 
.D m in,o,i- I n the case of D^, oj G C\{0}, one analogously obtains 

Ao G a p (D ul + B) if and only if — Ao G cr p (Dzj+ B), Ao G C\{0}, ui G C\{0}, 

(5.26) 

since 

(o 2)* Gt - b (o 2)=-^ s - ^ 

More precisely (cf. Lemma 5.6 below), the point spectrum of a(D + B) and hence 
that of «G T , a /p 2 : or equivalently, that of the associated quadratic pencil L(z) = 
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z 2 I-h + ziap~ 2 — T*T in leading order (independently of the boundary conditions 
used at x = 0, 1), consists of two infinite sequences {X±j(D + B)}j e fi C C, such 
that 



\ ±>j (D + B) = ±jn 



J 

Jo 



dx p{x) 



+ o(\j\). 



(5.28) 



Moreover, higher-order expansions of the type 



M 



X±,j(D + B) = ±jc_!+ Vc ro (±j)- ro + (|i| 

7 — < 4 



-M\ 



M e N, 



m=0 



C-l 



1 -i-l 



(5.29) 



have been studied under a variety of additional smoothness assumptions on the 
coefficients p and a, for instance, in [13], [14], [17], [23], [50], [78], [86], [87], [89], [90], 
[92], [110], and the references cited therein. Multiplicity questions of eigenvalues 
were discussed in [72] and [78] , and various eigenvalue inequalities were derived in 
[72]. 

Since the leading-order asymptotics for the non-self-adjoint operator D + B is 
of importance later in this section and we were unable to locate the result in the 
literature under our general assumptions on the coefficients a and p, we now present 
a proof of (5.28): 

Lemma 5.6. Assume Hypothesis 5.1. Then the leading- order asymptotics (5.28) 
holds. Consequently, the identical leading-order eigenvalue asymptotics applies to 
iGT.a/p 2 an d the associated quadratic pencil L(z) = z 2 I + ziap~ 2 — T*T, z £ 



C. Here D and T represent any of fmax,min,o,i,w and T, 
respectively. 



max, min 



o,i,u,, w e c\{0} ; 



Proof. First we note that by the results of Appendix A (cf. (A.3), (A. 10), (A.26), 
(A. 29), and (A. 30)) it suffices to focus on the nonnegative self-adjoint operator D 2 
(rather than on the Dirac-type operator D) and compare it to (D + B) 2 . Since 
D 2 = T*T 8 TT\ and T*T and TT* arc nonnegative self-adjoint Sturm-Liouvillc 
operators with appropriate boundary conditions at x — 0, 1 described in Section 3, 
the leading-order asymptotics of these operators is well-known to be independent 
of the boundary conditions used in this manuscript and given by (see, e.g., [Ill, 
Theorem 4.3.1]] 

-2 



Xj(T*T),Xj(TT*) = j 2 n 2 



dx p(x) 



Next, one chooses p > sufficiently large such that 



(D 2 + pI 2 )- l / 2 [BD + DB + B 2 ](D 2 + ph)- 1 ' 2 



(5.30) 



< 1, (5.31) 



B(L 2 ([0,l];p 2 dx) 2 ) 

where we recall our simplifying notation I2 — lL 2 ({o,i]-p 2 dx) 2 , an d observe that 



(D 2 + pI 2 )-^l 2 \BD + DB + B 2 }(D 2 + ph)- 1 ' 2 
= (D 2 + ph)- 1/2 [BD + B 2 }(D 2 + phT 1/2 
+ D(D 2 + pI 2 )- 1 / 2 B(D 2 + pI 2 )- 1 / 2 . 



(5.32) 
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Then the identity 
((£> + Bf + 



(D 2 + iih 



-1/2 



x / , + ( + /;/;;, ' • HI) + VB+ IS-' J) - + fj,l 2 ) ' • •: I) ' 
= (D 2 +l ,I 2 )-^ 2 [I 2 +C](D 2 +pI 2 )-^ 2 , 
where we introduced the Hilbert-Schmidt operator C in L 2 ([0, 1]; p 2 dx) 2 



(5.33) 
(5.34) 



C = 



h + (D 2 + nI 2 )-i/*[BD + DB + B 2 }(D 2 + iih)- 1 ' 1 



-h 



h + {D 2 + iiI 2 )-V*[BD + DB + B 2 ](D 2 + nI 2 )~V 2 



(5.35) 



x (D 2 + tiI 2 )-i/ 2 [BD + DB + B 2 ](D 2 + /iJ 2 )-V2 g £ 2 (L 2 ([0, 1]; p 2 dx) 2 ) 
recalling (L> 2 + /x/ 2 )" 1/2 G B 2 (£ 2 ([0, l];p 2 cfa) 2 ) (cf. (5.9)). Since clearly, 



! + (D 2 + iiI 2 )-^l 2 \BD + DB + B 2 ](D 2 + fih)- 1 / 2 

E B(L 2 ([0,l];p 2 dx) 2 

Theorem V.11.3 in [43] applies and hence 

lim Xj((D + B) 2 + pI 2 )/X 3 (D 2 + fil 2 ) = 1, 



(5.36) 



(5.37) 



and thus also lim.,-^ Xj((D + B) 2 )/X 3 (D 2 ) = 1. Together with (5.30) and the 
results in Appendix A relating the spectra of D and T*T, TT*, mentioned at the 
beginning of this proof, yield the leading-order asymptotics (5.28). □ 

Next, we recall a useful result due to Livsic [65] on the trace of the imaginary part 
of certain non-self-adjoint operators. But first we mention some preparations: Let 
F e Boo('H), then F has a complete system of eigenvectors and generalized eigen- 
vectors in % if the smallest linear subspace (i.e., the linear span) of all eigenvectors 
and generalized eigenvectors of F is dense in H. 

We continue with Schur's lemma: 

Lemma 5.7 ([42], Lemma II.3.3, [43], Theorem V.2.1). 

Let F G BoolTi) and denote by £(F) the smallest closed linear subspace of H (i.e., 
the closed linear span) containing all eigenvectors and generalized eigenvectors of F 
corresponding to all nonzero eigenvalues of F (i.e., to tr(F)\{0}). Then there exists 
an orthonormal basis {<Pj}j£j of £(F) (with J C N an appropriate index set), such 
that with respect to the basis {<Pj}jeJ, F is an upper-triangular operator satisfying 



with 



k=l 

Fjj = ^j(F), jeJ, 

a(F)\{0} = {X 3 (F)} 3eJ , \J\- 



(5.38) 



]T m(A, F). (5.39) 

Ae<r(F)\{0} 

Here |J| abbreviates the cardinality of J. The orthonormal basis {fj}je.J 01 
£(F) is obtained from building Jordan blocks associated with chains 

= *j(F)i>i,j, 



THE DAMPED STRING PROBLEM 



25 



Alp m ,j = Xj(F)lj) m j Or Alp m j = \j(F)l[) m ,j +lf) m -l,j, (5.40) 

m = 2,...,m a (X j {F),F), j e J, 

followed by the Gram-Schmidt orthogonalization procedure of the system ip m ,j, 
m = 1, . . . , m a (Xj (F), F), j G J. 

Next, we briefly consider bounded operators G G B(H) with nonnegative imag- 
inary parts, Im(G) > 0. Then (as shown, e.g., in [42, p. 136] in the case fio = 0), 
/ G ker(G — Holy) f° r some /i G M implies (/, (G — HoIn)f)n = 0, hence, 
(/, (G* - H I H )f)H = 0, and thus, < (f,lm(g)f) n = 0. Consequently, Im(G)/ = 
(cf. (A. 8)) and hence / G kcr(G* — HqI-h)- The symmetry of this argument with 
respect to G and G* yields 

kcr(G - nol-u) = kcr(G* - (i I n ), ^ G R, (5.41) 

and hence also 

H = kcr(G - nal n ) ® ran(G - ^I H ) (5.42) 



(rather than the standard % = ker(G* — plqI-h) © ran(G — HqI-h))- Relation (5.42) 
is interesting as it implies the following fact: 

If ker(G — Holy.) ^ {0}, then G has no generalized (resp., associated) 
eigenvector corresponding to the eigenvalue Ho EM.. 

Indeed, assuming that for some /o G H, 

(G - /i /«) 2 /o = but (G - MoWo + (5.44) 

then yields 

(G - HoI H )f G [kcr(G - n Q I u ) D ran(G - Mo J„)] (5.45) 

and hence contradicts (5.42). 

The following trace formula (5.51), the centerpiece of the next theorem, is a 
well-known result due to Livsic [65] : 

Theorem 5.8 ([42], Lemma VII.8.2, [43], Lemma II.4.1). 

Assume that G G S (X) ('H) and Im(G) > 0. In addition, denote by £{G) the smallest 
closed linear subspace (i.e., the closed linear span) of all eigenvectors and generalized 
eigenvectors of G corresponding to the nonzero eigenvalues of G . 

(i) Then £(G) C ran(G). Moreover, 

G has a complete system of eigenvectors and generalized eigenvectors (5.46) 
if and only if 

£{G) = ran(G), (5.47) 

which in turn is equivalent to 

8(G) 1 - = ker(G). (5.48) 

(ii) Suppose in addition that Im(G) G B\(%). Then 

^Im(Ai(G))<tr w (Im(G)). (5.49) 

jeJ 

Here {Xj(G)}j^j, J C N an appropriate index set, denotes the eigenvalues of G 
ordered with respect to decreasing magnitude, 

|A,+i(G)| < |A,(G)|, j,j + leJ, (5.50) 
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counting algebraic multiplicities (with \J\ = X)ago-(g)\{o} m (X, G)). In addition, 

^Im(A j (G))=tr«(Im(G)) (5.51) 

jeJ 

if and only if 

G has a complete system of root vectors (i.e., eigenvectors and 

generalized eigenvectors). (5.52) 

For convenience of the reader we present the short argument for the inequality 
(5.49) and indicate the equivalence of (5.51) and (5.52): Choosing an orthonormal 
Schur basis {<j>j} jeJ for G in £(G) (cf. (5.38), (5.39)), 

3 

= Gkj'Pk, j G J, (5.53) 
k=i 

where 

Gj,j = {(f>j,G4>j) n = Xj(G), G fc j = (<f> k , G4>j) H , j, k G J, (5.54) 

and an orthonormal basis {xk}k£K of 6(G)- 1 , with K C N an appropriate index 
set, one obtains 

tr„(Im(G)) = ^(^,Im(G)^)„ + £ ( Xk , Im(G) Xfc )w 

= E Im (A 3 (G)) + ^(xfe,Im(G) Xfe )w, (5.55) 

using 

(^,Im(G)0,)„ = Im((^,G0,)„) = Im(Aj(G)), j G J. (5.56) 
Since by hypothesis Im(G) > 0, (5.55) proves (5.49). Moreover, (5.51) holds if and 
only if lm(G)\g(Q}± = 0. The latter can be shown to be equivalent to S(G) ± = 
ker(G), which in turn is equivalent to £(G) = ran(G) by (5.42). 

Since Theorem 5.8 (i) requires completeness of the system of eigenvectors and 
generalized eigenvectors of G, we next recall a sufficient criterion for completeness 
convenient for our subsequent purpose: 

Theorem 5.9 ([42], Theorem XVII.5.1, [68], Theorem 1.4.3). 
Let S be a self-adjoint operator in H with purely discrete spectrum, or equivalently, 
satisfying (S — zly)^ 1 G Boo(Tl) for some (and hence for all) z G p(S). Moreover, 
let T be an S-compact operator in %, that is, T(S — zlu)^ 1 G B 00 ('H) for some 
(and hence for all) z G p(S). In addition, denoting by {Xj(S)}j e ,j, J C N, the 
eigenvalues of S ordered with respect to increasing magnitude, \Xj(S)\ < \Xj+i(S)\, 
j G J, counting algebraic multiplicities, assume that for some p > 1, 

\^(S)\~ p < oo, (5.57) 

jeJ' 

where J' C J represents the index corresponding to all nonzero eigenvalues of S . 
Then a(S + T) consists of only discrete eigenvalues in the sense that 

(S + T-zIu)- 1 eBooiH), ze P (S + T), (5.58) 

and each element Xk(S + T), k G K , K C N, of a(S + T) is an isolated point of 
a(S + T), each Xk(S + T) has finite algebraic multiplicity (i.e., the range of the 
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Riesz projection associated with Xk(S + T) is finite- dimensional) , and the system 
of eigenvectors and generalized eigenvectors of S + T is complete in H. 

In particular, the system of eigenvectors and generalized eigenvectors of (S + T — 
zl-ji) -1 associated with all (necessarily nonzero, cf. (5.58)) eigenvalues of (S + T — 
zlu)^ 1 is complete in H. We also note that the relative compactness assumption 
on T with respect to S yields that 

\\T(S - zI n )~ 1] \ B{;H) < 1 for z G p(S), \z\ sufficiently large, (5.59) 

and hence 

(S + T- zlu)- 1 = [I H + T(S - z)- 1 } ~\s - zlu)- 1 G B^H), (5.60) 

for z G p{S), \z\ sufficiently large, is well-defined (cf. [104, p. 200 and Theorem 
9.7]), implying (5.58). 

Remark 5.10. The proof of Theorem 5.9 in [42, Theorem XVII.5.1] and [68, Theo- 
rem 1.4.3] relies on a completeness result of Keldysh (reproduced in English in [42, 
Theorem IX. 4.1] and [43, Theorem V.8.1]) in the context of certain multiplicative 
perturbations of compact operators. In addition, we note that Theorem 5.9 extends 
to a normal operator S whose spectrum lies on a finite number of rays starting at 
z = 0, moreover, it suffices to take p > (cf. [68, Theorem 1.4.3]). We also re- 
mark that Theorem XVII.5.1 in [42] requires S^ 1 G B 00(H), TS' 1 G B 00(H), and 
(in+TS^y 1 G B(H), but if G a(S), the simple replacement of S by S- z I n , 
for appropriate ^ z G R, permits one to remove the restriction of bounded 
invertibility of S (cf. (5.58)). 

Combining Theorems 5.3, 5.4, 5.9 and Theorem 5.8 (i) then yields the following 
infinite sequence of trace formulas for the damped string equation: 

Theorem 5.11. Assume Hypothesis 5.1 and denote by {Xj(D+B)}j e j, JCZ, the 
eigenvalues of D + B ordered with respect to increasing magnitude, \Xj(D + B)\ < 
\Xj+i(D + B)\, j G J, counting algebraic multiplicities. In addition, assume that 
G p(D + B). Here (and below) T, D, andt 2n , n G No, stand for one o/T m i ni o,i, w , 
-Dmin,o,i.w; and tmin,o,i,aj;2n; n G No, respectively, and we only consider the case 
ui G C\{0, 1} (cf. (3.20)). Then the following infinite sequence of trace formulas 
hold: 

^ Im(A J -(D + gr+ 1 ) j-t 2n , m = 2n, 

g \X j{ D + B)^ ={ , TO = 2n + 1 , ^N„. (5.61) 
Explicitly, one obtains for m = 0, 1 in (5.61), 

\\-(D + B)\ 2 ~ ~ trL2 ([o,i];p 2 dx)(WP )(T T) ) 



- dxx(l — x)a(x), 

- Jq dx xa(x), 

- dx (1 — x)a(x), 

k jjE^ Jo ^ [(1 - \u\ 2 )x + l]a(x), uj G C\{0, 1}, 



T 




T 


= T Q , 


T 


= T U 


T 


— T 



(5.62) 
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v MA j (Q + fl))R C (A j (fl + 5)) 

£ |A,(^ + B)| 4 -°" ^ 

Proof. By (5.9) and (5.10), the hypotheses of Theorem 5.9 are satisfied, identifying 
S and P with D and £?, respectively. In particular, a{D + B) consists of only 
discrete eigenvalues, {\j(D + B)}j e j, J C Z, and the corresponding system of 
eigenvectors and generalized eigenvectors of D + B is complete in L 2 ([0, 1]; p 2 dx) 2 . 

However, the crucial hypothesis Im(G) > in Theorem 5.8 is not necessarily 
implied by Hypothesis 5.1. To remedy this fact one can proceed as follows: The 
operators D + B and D + B + zl 2 , z E C, have the same set of eigenvectors 
and generalized eigenvectors which are complete in L 2 ([0, l];p 2 dx) 2 . As in (5.12) 
one concludes that for \(\ < e , < £o sufficiently small, using that B is purely 
imaginary, Re(B) = 0, cf. (5.4), 

ImliD + B-iz + Oh)- 1 ] 

= (D + B-(z + Ohy^lmiz + C) - Im(S)] [(D + B - (z + Oh)' 1 ] * > 0, 

Im(z) > £ + ||-B|| B(£ 2 ([0)1];p 2 da . ) 2 ) . (5.64) 

Indeed, recalling the Hilbert-Schmidt property (5.9) of the resolvent of D + B, one 
infers that (5.64) is well-defined since 

(D + B - (z + Oh)- 1 =(D-(z + Oh)- 1 

x [h+BiD-iz + OhrT 1 eB(L 2 ([0,l];p 2 dx) 2 ), 

lm(z) > £ + ||-B|| B ( L 2 ([0il];p 2 da . ) 2 ) , (5.65) 

and since 

\\B(D- (z + Oh) 1 || B ( i 2 ([0i i ];<:) 2 dx) 2 ) 

< ll-B||f5(L2([04];p 2 d a; )2)[Ini(z) - £ ] _1 < 1, 

as Im(z) > £ + \\B\\ B{n) . 

At this point one can apply Theorem 5.8 (ii) (especially, (5.51), (5.52)) to obtain 

tr L 2 ([0;1] .„2 dx) 2 (Im ((£> + B - (z + Oh)' 1 )) 

= J2 Im (( x i( D + B ) -z-Cy 1 ) < (5-67) 

|C| < £o, Im(z) > £ + ||-B|| B(W) . 

Since (5.67) exhibits no analyticity with respect to z we cannot simply continue it 
to z = 0. As a result we need to proceed along a different route: First we note that 
e p(D + B) and the asymptotic behavior (5.28) of the eigenvalues of D + B yield 

£|A,-(i? + B)|- 2 <oo. (5.68) 
jeJ 

Recalling (5.24), this also implies 

SIW + *)-)|=E ^i^ <~- ("») 

je.J jeJ 1 Jy " 

Thus, choosing a fixed Im(zi) > £ + ||B|| B ( W ), and using the Hilbcrt-Schmidt 
property (5.9) of the resolvent of D + B once more, one computes with the help of 
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(5.67) that 

^L 2 ([0.1];p 2 dx) 2 (lm((£> + J B-C/ 2 )- 1 )) 

= tr L 2 {m . p 2 dx) 2 ( Im ((£> + B - ( Zl + Oh)' 1 )) 
+ trL2 ([0>1] .„2 dx) 2 (Im ((D + B - (h)' 1 )) 

- tr L 2 ([0il] . p 2 da:) 2 (Im ((£> + B - (« a + O-M -1 )) 
- tr L 2 ([oa];p 2 d;c) 2 (Im ((£> + B - (z x + O-fc) -1 )) 

+ tr z ,2 ([0;1] . (9 2 (jx) 2 (Im ({D + B - C^) -1 ) 

-TmdD + B-izi+Qh)- 1 )) 
= tr L 2 m] . p 2 dx) 2 (lm((D + B- ( Zl + Oh)' 1 )) 

+ Im (tr L2([04];p 2 dx)2 ((£> + B - (hV 1 — (D + B — (z\ + Qh)' 1 )) 

= ^Im((A,p + S)-z 1 -C)- 1 ) 
jeJ 

- Im L ]T(A, (D + B) — zi — C)- 1 (Xj (D + B)-Cy 1 ) 

= J2lm((X j (D + B)-C)- 1 ), ICI < ^0- (5.70) 

Thus, one obtains for ( G (— £0, £o)\{0}, < £0 sufficiently small such that 
(-eo.eo) C p(D + B), 

tr i 2 ([0 , 1];p 2 dx) 2 (im ((D + B — (h)- 1 )) = £lm ((Xj(D + B) - C)" 1 ) 

= ^Im (A^D + £?) _1 (l - CA,p + B)- 1 )- 1 ) 
jeJ 

= E Im ( E A ^ + B ) _1 (cw + s)- 1 )") 

= ^^Im(A j ( J D + B)-"- 1 )C" 

= E (EMW+sr^V 

v / v Im(A J (D + i?)»+ 1 ) \ 

= E *2nC 2 ", (5-71) 
nSNo 

proving (5.61) subject to the interchange of the sums over j and n. To justify this 
interchange it suffices to prove the absolute convergence of X^jej S n eN ^ m 
-S) - ™ 1 )C™- F° r this purpose one uses 

n 

Im(z" +1 )=Im(z)^z l I , '- f , zeC, (5.72) 
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and estimates 

n 

|lm(z™ +1 )| < |Im(z)|^|z| £ |zr i - £ = (n + l)|Im(z)||zr, z G C. (5.73) 

1=0 

By (5.69) one estimates, without loss of generality, assuming that \ J\ = oo (other- 
wise, there is nothing to prove) and for simplicity, using J = Z\{0}, 

N N 

]T ^llm^OD + i?)-"- 1 )^ 

j=-N n=0 

j = -Nn=0 1 jK 

\j=-N j=M(s)/ n=0 1 • JV 71 

ICI<£, 

for some C(e) > 0, < M(e) < iV, and some C > independent of TV G N, proving 
the interchangeability of the sums over j and n. Here we used that \Xj(D + B) \ — > 
oo as \j\ — > oo and, since by hypothesis G p(D + B), \Xj(D + B)\ > e, and 
\\j(D + B)\~ n < £-", j G J n [M(e), oo), for sufficiently small e > 0. 

Taking m = 0, 1 in (5.61) and applying (5.16)-(5.19) then yields (5.62) and 
(5.63). □ 

Remark 5.12. (i) Trace formulas, quite different from the ones we derived in Theo- 
rem 5.11, involving certain regularized sums of all eigenvalues were derived in [13], 
[14], [17], and [110]. For instance, in the case of Dirichlet boundary conditions at 
x = 0, 1, and for p(x) = 1, the trace formula 

^[A_ j{D min + B) + A +J (A„i„ + B) - 2c ] - l -[a{Q) + a(l)] + c , 

(5.75) 

co = / dxa(x), 
2 Jo 

was derived in [13]. 

(ii) We emphasize that dissipativity of B, more generally, sign-defmitcness of a, is 
not assumed in this section. 

Next we briefly turn to the question whether the system of (algebraic) eigen- 
vectors associated with D + B as in Theorem 5.11 constitutes a Riesz basis. One 
recalls that a system {/„}„ g n in % represents a Riesz basis in H if it is equivalent 
to an orthonormal basis {e„}„ e N in %, that is, if there exists an operator V G B(H) 
with V~ x G B(H) and Vf n — e„, n G N. Equivalently, for any h G W, there exists 
a sequence {c n (/i)}„ e N C C, such that 

h = ^2c n (h)f n (5.76) 

riGN 

converges unconditionally in T-L (i.e., it remains convergent to the same element in 
TL under any permutation of the terms c n (h)f n , n G N, in (5.76)). 
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More generally, a system {g n }nen in H represents a Riesz basis with parentheses 
in H if there exists a sequence of strictly increasing positive integers {mi}ien C 
N, 1 < mi < me + i, I G N, such that for any h G H, there exists a sequence 
{(i„(/i)}„ e N C C, implying the unconditional convergence of the following series 
over ^ G N in H, 

(m e+ i \ 
]T d n (h)g n = ]T P t h, (5.77) 

where m = and the projections Pi, i G N , are given by 

m e+1 

P t h= d n(h)g n , he%, P e P m = P m Pe = S e , m Pm, e,meN (5.78) 

n— + 1 

(with mi, I e No, independent of h G 14.) . Equivalently, the sequence of subspaces 
{ran(P£}^ eNo with P t satisfying PiP m = P m Pi = Si, m P m , l,m € N , is called an 
unconditional (or Riesz) basis of subspaces in H. 

The following fundamental abstract result on the existence of a Riesz basis with 
parentheses due to Katsnclson [62], [63], Markus [67], and Markus and Matsaev [69], 
[70] (proved under varying generality) will subsequently be applied to the operator 
D + B in L 2 ([0,\];p 2 dx) 2 : 

Theorem 5.13. Let N be a normal operator in H with compact resolvent and 
spectrum lying on a finite number of rays lZ q = {z G C | arg(z) = 8 q G [0, 27r)}, 
q = 1, . . . , r for some r G N. In addition, suppose that for some p G [0, 1), 

liminf^-^^iV) < oo, (5.79) 

tfoo 

where #(t,N), t > 0, denotes the sum of the multiplicities of all eigenvalues of N 
in the open disk D(0;t) C C of radius t > centered at the origin z = 0. Assume 
that R is a densely defined operator in H satisfying for some c > 0, 

dom(iV) C dom(i?), \\Rf\\ B(H) < c\\Nf\\ P B{n) \\f\\^ p , f G dom(iV) (5.80) 

(i.e., R is p- subordinate to N, cf. [68, Sect. 5] and hence R(N — zl-u)^ 1 G Boo(7i), 
z G p(N)). Then N + R defined on dom(N + R) = dom(A^) is a densely defined 
closed operator with compact resolvent and the system of root vectors of N + R 
forms a Riesz basis with parentheses in H. 

We note that since N is assumed to be normal, the geometric and algebraic 
multiplicty of all eigenvalues of N coincide in Theorem 5.13. 

For a detailed discussion of Theorem 5.13 we refer to Markus [68, p. 27-37, 
Theorem 6.12]. For additional results in connection with Theorem 5.13 we also 
refer to [1], [2], [3], [4], [5, Ch. 5], [10], [18], [29, Ch. XIX], [48], [49], [61, Sect. 
V.4.5], [64], [82], [83], [84], [85], [95], [102], [105], [106], [107], [108], [109], [112], and 
the references cited therein. 

Theorem 2.3 (especially, Theorem 2.4) permits one to shift spectral considera- 
tions of Ga,r in Ha®Hi or that of the quadratic pencil L(-) in (2.57), to that of the 
simpler Dirac-type operator Q + S in H (B H. In particular, in the concrete context 
of this section, B is a bounded diagonal perturbation of the supersymmetric self- 
adjoint Dirac-type operator D in L 2 ([0, 1]; p 2 dx) 2 and an application of Theorem 
5.13 yields the following result: 
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Theorem 5.14. Assume Hypothesis 5.1 and introduce the super symmetric Dirac- 
type operator D, the non- self- adjoint diagonal perturbation B, and the operator 
D + B in L 2 ([0,l];p 2 dx) 2 by (5.3), (5.4), and (5.5), respectively. Then D + B 
is closed and densely defined on dom(D + B) = dom(D) with compact resolvent. 
Moreover, the system of root vectors of D + B, or equivalently, that for Gt,o/p 2 
(resp., the quadratic pencil N(z) = z 2 I + ziB — T*T , dom(A(z) = dom(T*T), z E 
C), forms a Riesz basis with parentheses in L 2 ([0, 1]; p 2 dx) 2 (resp. L 2 ([0, 1]; p 2 dx)). 
More precisely, each of the two sequences {\±j(D + B)}j e iq in (5.28) decomposes 
into finite clusters of eigenvalues A± t t tending to ±oo as I — > ±oo such that the 
associated Riesz projections P(A±j, D + B), I E No, satisfy 

J2 P{A—,t, D + B)+J2 p ( A +,e, D + B)=I. (5.81) 

£en„ ien a 

Here D and T stand for one of -D m ax,min,o,i,w and T maX)m i n) o,i lU ;, respectively. 

Proof. That D + B defined on dom(D + B) = dom(£>) is a densely defined closed 
operator with compact resolvent is clear since B E S(i 2 ([0, 1]; p 2 dx) 2 ) and the 
resolvent of D is compact in L 2 ([0, 1]; p 2 dx) 2 . In order to apply Theorem 5.13 we 
first note that by (5.25) all eigenvalues of D + B lie in a strip symmetric with respect 
to the real axis and that by (5.25) and (5.28) the eigenvalues of D + B split up 
into two sequences converging to ±oo. As a result of the proof of Theorem 5.13 as 
described in detail in [68, p. 27-37], one can now apply Theorem 5.13 in the special 
case where N = D is self-adjoint, r = 2, R = B is bounded, and one chooses p = 
since #(i, D) = 0(t) as 1 1 oo, given the asymptotics (5.28). □ 

We conclude with the following observations: 

Remark 5.15. It appears that Hypotheses 5.1 is the most general set of assumptions 
under which the Riesz basis property with parentheses has been proven for D+B, or 
equivalently, G T a j p i (resp., the quadratic pencil N(z) = z 2 I + ziB — T*T, z E C). 
Under additional smoothness assumptions on a and p, and typically for separated 
boundary conditions in T*T, TT* , and/or in the case of damping at the end points 
x = 0,1, the Riesz basis property (without parentheses) for Gt,o/p 2 has been 
established in [23], [24], [47], [86], [87], [88], [90], [91], [92]. In this context we recall 
that the existence of a Riesz basis (without parentheses) for D + B, or equivalently, 
Gr.a/p 2 , is equivalent to both operators being unbounded spectral operators in 
the sense of Dunford (cf. [29, Chs. XVIII, XIX]). In the case of the boundary 
conditions T max ^ m i n ,o,i,w, w G C\{±1}, one can expect the Riesz basis property 
(without parentheses) to hold under our general assumptions in Hypothesis 5.1. 
For a particular class of one-dimensional Dirac-type operators the delicate question 
of unconditional convergence of spectral Riesz expansions has recently been treated 
in great detail by Djakov and Mityagin [26], [27], [28], [71], with earlier contributions 
in [21], [52], [73], [98], [99]. 

We also note that scalar spectral operators were discussed in connection with 
(generalizations of) damped wave equations by Sandefur [82]. 

Appendix A. Supersymmetpjc Dirac-Type Operators in a Nutshell 

In this appendix we briefly summarize some results on supersymmetric Dirac- 
type operators and commutation methods due to [25], [40], [96], and [97, Ch. 5] 
(see also [51]). 
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The standing assumption in this appendix will be the following: 
Hypothesis A.l. Let Hj, j = 1,2, be separable complex Hilbert spaces and 

A : dom(A) QHi^H 2 (A.l) 
be a densely defined closed linear operator. 

We define the self-adjoint Dirac-type operator in Hi © H 2 by 
'0 A* 



A 



dom(Q) = dom(A) © dom(A*). (A.2) 



Operators of the type Q play a role in supersymmetric quantum mechanics (see, 
e.g., the extensive list of references in [12]). Then, 



AA\ 

and for notational purposes we also introduce 

Hi = A* A in Hi, H 2 = AA* in H 2 . (A.4) 

In the following, we also need the polar decomposition of A and A*, that is, the 
representations 

A = V A \A\ = \A*\V A , A* = V A ,\A*\ = \A\V A », (A.5) 

where 

|A| - (A*A) 1/2 = h\ /2 , \A*\ = (AA*) 1/2 = Hi 12 . (A.6) 

Here V A is a partial isometry with initial set ran(|A|) and final set ran(^4) and V A * 
is a partial isometry with initial set ran(|A*|) and final set ran(A*). In particular, 



J A(A'A)-V* = (AA*)~V*A on (ker(A))\ 
VA = (0 on ker(A). (A ' ?) 

Next, we collect some properties relating Hi and H 2 : 

Theorem A.2 ([25]). Assume Hypothesis A.l and let <p be a bounded Borel mea- 
surable function. 

(i) One has 

kei(A) = ker(fTi) = (ran(A*)) ± , kcr(A*) = kcr(ff 2 ) = (ran(A)) ± , (A.8) 

V A H n J 2 = H r 2 l/2 V A , n e N, V^(#i) = ^JVa- (A.9) 

(ii) Hi and H 2 are essentially isospectral, that is, 

<t(Hi)\{0} = o-(H 2 )\{0}, (A.10) 

in fact, 

A* All-Hi — Pkcr(A)} i s unitarily equivalent to AA*[If{ 2 — Pker(A*)]- (A-ll) 
In addition, 

f e dom(JTi) and H x f = A 2 /, A ^ 0, 

implies Af e dom(# 2 ) and H 2 (Af) = \ 2 {Af), (A.12) 
<7 G dom(i? 2 ) and H 2 g ~ [i 2 g, /i^O, 

imp/ies G dom(#i) and H X {A* f) = n 2 {A*g), (A.13) 
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with multiplicities of eigenvalues preserved. 
(Hi) One has for z G p{H\) R p(H 2 ), 

I H2 + z(H 2 - zlu,)- 1 D A(Hi - zI Hi y x A\ (A.14) 

I Ul + z(R x - zlu,)- 1 D A*(H 2 - zln^A, (A.15) 

and 

A*<f>(H 2 ) D <j>{Hi)A*, Aij>{Hi) 2 4>(H 2 )A, (A.16) 

Va.^H^D^HJVa., VA^Hi) D <f>(H 2 )V A . (A.17) 

As noted by E. Nelson (unpublished), Theorem A. 2 follows from the spectral 
theorem and the elementary identities, 

Q = V Q \Q\ = \Q\V Q , (A.18) 

ker(Q) = ker(|Q|) = kcr(Q 2 ) = (ran(Q))- 1 , (A.19) 
I nM + z(Q 2 - zlu^y 1 = Q 2 (Q 2 - zln^u,)- 1 D Q(Q 2 - zI^n^Q, 

z e p(Q 2 ), (A.20) 

QHQ 2 ) 2 4>{Q 2 )Q, (A.2i) 

where 

(Va)*\_{0 V A .\ (A22) 



Ja ) \V A 

In particular, 

'1\ 

hcer(A*) 

and we also recall that 



kcr(Q) = ker(A) © ker(A*), P kcr(Q) = ( J k ^ A > p (A.23) 



<T 3 Qa 3 = -Q, c-3 = ^ , (A.24) 

that is, Q and — Q are unitarily equivalent. 

Finally, we note the following relationships between Q and Hj, j = 1,2: 

Theorem A. 3 ([12], [96]). Assume Hypothesis A.l. 
(i) Introducing the unitary operator U on (ker(Q)) ± by 

u = 2 -i/2 (Itu [Va)*\ on(ker(Q)) ± ) (A . 25) 



one infers that 



(ii) One has 



H 1 ' 2 



UQU- 1 = [ i 1/2 ) on (ker(Q))^. (A.26) 



(Q C/WlffiwJ -U^i-C 2 ^)- 1 C^-C 2 ^)" 1 J' (A.27) 



C 2 ep(i?i)n# 2 ). 
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(Hi) In addition, 

(£) edo^Q) (£)=,(£), ,^0, (A2g) 

implies fj <E dom(Hj) and Hjfj = rj 2 fj, j = 1,2. 

Conversely, 

f G dom(Jfi) and #i/ = A 2 /, A ^ 0, 

( f \ ( f \ ( f \ ( A ' 29 ) 

implies I . f \ 6 dom(Q) and Q . _{ . , = A J 



Similarly, 

f" ,:„n,^ and ^" J ^ ^ 
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